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Abstract 


A prediction method has been written and incorporated into a three- 
dimensional Navier-Stokes code (PAB3D) for the calculation of nozzle 
internal performance. The quantities calculated are discharge coeffi- 
cient; normal, side, and axial thrust ratios; rolling, pitching, and yawing 
moments; and effective pitch and yaw vector angles. Four different case 
studies are presented to confirm the applicability of the methodology. In- 
ternal and, in most situations, external flow-field regions are required to 
be modeled. The computed nozzle discharge coefficient matches both the 
level and the trend of the experimental data within quoted experimen- 
tal data accuracy (0.5 percent). Moment and force ratios are generally 
within 1 to 2 percent of the absolute level of experimental data, with the 
trends of the data matched accurately. 


Introduction 

Highly maneuverable aircraft operate over power 
settings and Mach numbers that require a propulsion 
system with variable geometry for obtaining efficient 
performance at different throttle settings. Under- 
standing the effects of various nozzle geometries on 
internal flow and on the surrounding boattail-nozzle 
region is vital for designing efficient afterbodies for 
these aircraft. The development and utilization of 
advanced computational methods plays a vital role 
in developing this understanding. Presently, there 
are several ongoing research activities at the Langley 
Research Center directed at establishing an experi- 
mental data base for new nozzle concepts. Subse- 
quent improvements to the computational methods 
are guided by these data. 

A nozzle internal performance module has been 
written, and the prediction capabilties of this mod- 
ule have been evaluated for a series of test cases. 
The module was incorporated into the Navier-Stokes 
solver PAB3D (ref. 1), which provides the flow-field 
solution for the prediction subroutine. The control 
volume concept was used for the calculation of the 
resultant body forces due to the fluid flow exiting 
the nozzle (ref. 2). Forces and moments are calcu- 
lated from the integration of the momentum fluxes 
through the volume faces. These integrations can be 
performed at intermediate steps throughout the so- 
lution procedure to provide a flow convergence quan- 
tity and to monitor the flow code’s computational 
residuals. 

This paper presents four nozzle case studies for 
evaluation of the performance prediction method- 
olgy. Each nozzle geometry has characteristics that 
require the calculation of different flows and perfor- 
mance quantities for each case. Comparisons of com- 
puted discharge coefficient, axial thrust ratio, resul- 


tant, vector angles, and pitching and yawing moments 
with experimental data are presented as appropriate 
for the following cases. 

Case 1: The first case study is from the series 
of Stratford choke nozzles (ref. 3) that serves as 
the calibration standard for many propulsion testing 
facilities. These geometries are convergent nozzles 
with a circular-arc internal closure contour through 
the throat. These nozzles have well-established mass- 
flow characteristics that provide a highly accurate 
standard for comparisons of performance predictions. 
A grid convergence study of discharge coefficient is 
shown for the nozzle with a throat area of 25.766 cm 2 
(3.992 in 2 ) operating at a pressure ratio of 2. 

Case 2: The capability of predicting discharge 
coefficient and axial forces for a type of convergent- 
divergent (C-D) nozzle is verified with the experi- 
mental data of Mason, Putnam, and Re (ref. 4). 
Their data were obtained on a model designed to 
address the performance characteristics due to para- 
metric changes of nozzle internal geometry. The noz- 
zle throat contour, that is, the radius of the circular- 
arc throat contour, was systematically changed in a 
series of nonaxisymmetric C-D nozzles with parallel 
sidewalls. Internal static pressures, thrust ratio, and 
discharge coefficient data for various nozzle operat- 
ing pressure ratios were measured. Low- and high- 
expansion-ratio geometry nozzles were investigated. 
For this study, comparisons of prediction with ex- 
perimentally determined nozzle discharge coefficient, 
thrust ratio, and centerline divergent flap static pres- 
sures are presented for the low-expansion-ratio noz- 
zles only. 

Case 3: The third geometry provided data to 
evaluate the prediction of side forces coupled with 
axial forces. The translating sidewall yaw vector 
concept tested by Mason and Berrier (ref. 5) was one 


of several yaw vector concepts on which discharge 
coefficient, thrust ratio, resultant thrust ratio, and 
effective pitch and yaw vector angles were measured. 
Comparisons of prediction with an experimentally 
determined effective yaw vector angle are presented. 

Case 4: The axisymmetric nozzle with multiaxis 
thrust vectoring that was tested experimentally by 
Carson and Capone (ref. G) provided data for evalu- 
ation of the prediction of pitching moments coupled 
with yawing moments. Five combinations of geomet- 
ric* pitch and yaw vector angles were tested in dry 
power and afterburning power settings for two dif- 
ferent lengths of divergent flaps. The nozzle expan- 
sion ratio was 1.35 for both settings, and thus the 
design nozzle pressure ratio w r as 5.01. A dry power 
setting would be a throat area permitting a paiticu- 
lar mass How indicative of aircraft cruise thrust lev- 
els. The data available include discharge coefficient 
and forces, moments, and internal static pressures at 
several angular stations around the nozzle. Only per- 
formance data from the short-flap, dry power nozzle 
are presented. 

Symbols and Abbreviations 

4 nozzle cross-sectional area, cm 2 

A r nozzle exit area, cm" 

,4/ nozzle throat area, cm 2 

A* nozzle cross-sectional area at sonic* 

conditions, cm 2 

AR nozzle aspect ratio, ratio of width to 
height at exit 

a local speed of sound, in/sec 

C d discharge coefficient. ^ 

e energy, J 

F gross thrust along body axis, N 

F total force vector, N 

F- ideal iscntropic gross thrust along body 

axis, N 

l length from nozzle connect station to 

nozzle exit station, cm 

[ s axial length from nozzle throat to nozzle 

exit, cm 

M total moment vector, N 

M (l iscntropic nozzle exit Mach number 

Mi local jet Mach number 


N P t 

NPR 

NPR„ 

n 

P 

Pt 

PtJ 

Poc 

R 

R 

TV 

T 

T t 

Tlj 

U 

u 

U V 

m 

w P 


Xt 

y + 


AA 

bv,p 

fiv.y 

6y 

e 

0 

v 
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free-stream Mach number 

number of grid points 

nozzle pressure ratio, 

design nozzle 1 pressure ratio 

surface normal vector 

static pressure, Pa 

total pressure. Pa 

jet total pressure, Pa 

free-stream static pressure. Pa 

gas constant (7 = 1.4), 287.3 J/kg-K 

moment arm vector, cm 

nozzle circular-arc throat radius, cm 

static temperature, K 

total temperature, K 

jet total temperature. K 

total velocity vector, in/sec 

velocity in streamwise direction, in /sec 

velocity at edge of boundary layer, m/sec 

ideal mass-flow rate, kg/sec 

calculated mass- flow rate, kg/sec 

axial distance measured from nozzle 
throat, positive downstream, cm 

axial distance from nozzle connect 
station to end of left sidewall, cm 

axial distance from nozzle connect 
station to nozzle throat station, cm 

nondimensional distance from w r all in a 
turbulent shear layer 

vertical ordinate, cm 

ratio of specific heats, 1.4 for air 

incremental cross-sectional area, cm 2 

pitch vector angle, deg 

yaw vector angle, deg 

effective yaw vector angle, deg 

nozzle expansion ratio, ^ 
nozzle convergence flap angle, deg 
kinematic viscosity, m 2 /sec 
density, kg/nV* 
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Experimental Data 

A sketch of a representative test model is pre- 
sented in figure 1. The Stratford nozzle geometry is 
shown in figure 2. Discharge coefficient data for the 
Stratford nozzle were obtained from the static thrust 
tables in the data reduction system used by the Lang- 
ley 16-Foot Transonic Tunnel. The other three nozzle 
geometries chosen for this ('valuation were tested in 
the static test facility of the 16-Foot, Tunnel. The 
two nozzles reported in references 4 and 5 were non- 
axisymmetric C-D type designs with circular-arc 
throat contours, straight divergent flaps, and flat 
nondiverging sidewalls; the relevant geometric pa- 
rameters are shown in figures 3 and 4. The nozzle 
reported in reference 6 (case 4) was an axisymmctric 
C-D type with a sharp throat corner, as shown in fig- 
ure r5 for the unvectored geometry. Figure 6 presents 
sketches of this nozzle at several thrust-vectoring an- 
gles. For each nozzle, the jet was simulated by high- 
pressure air exhausting out the exit of the nozzle into 
static air. The nozzle discharge coefficient was de- 
termined from the experimentally measured jet total 
temperature, jet total pressure, and mass-flow rate. 
Thrust ratio was determined from the measured bal- 
ance axial force that was nondimensionalized by the 
ideal thrust (determined from the measured nozzle 
mass flow). Multiaxis force measurements were de- 
termined from a six-component force balance. 

Computational Procedure 

Flow-Field Calculation 

A three-dimensional Nbivier- Stokes code PAB3D 
has been under development with specific applica- 
tions for predicting the effects of jet exhaust plume 
on three-dimensional (3-D) nozzle-afterbody config- 
urations. The thin-layer Navier-Stokes formulation 
(ref. 7) was modified to simulate jet mixing problems 
(ref. 1). The code allows for the partitioning of the 
flow-field domain into multi block grids and is capa- 
ble of using several numerical schemes to solve the 
governing equations and different turbulence models 
(ref. 8). 

A computational domain can be divided into 
zones consisting of longitudinal planes that are ar- 
ranged linearly in the streamwise direction. Each 
zone is subdivided into blocks that can have different 
boundary conditions applied at each of their six faces. 
The four lateral block faces either communicate with 
adjacent block faces or serve as outer zone boundaries 
or symmetry planes. The upstream and downstream 
faces are zone interfaces that are handled similar to 
the lateral faces. A user-written control file deter- 
mines the communication between blocks and the 


type of boundary condition to be used at each face. 
Blocks can be partitioned along the streamwise direc- 
tion to provide for a change in boundary conditions 
part way along a lateral block face. Different, grid 
topologies for neighboring blocks are permitted, with 
sonu' restrictions on grid matching at block bound- 
aries. This structure allows relatively complex con- 
figurations and flow conditions to be modeled more 
efficiently. Additionally, the code permits different 
numerical schemes to be applied selectively to each 
block, for example, the space-marching technique for 
supersonic flows and time-dependent calculations for 
subsonic and separated flows. 

Boundary Conditions 

The PAB3D cock' allows different, types of bound- 
ary conditions to be applied to any given block face. 
Solid walls arc treated as no-slip adiabatic surfaces, 
and the block interfaces are CO continuous. The 
boundary conditions used for the internal nozzle flow 
path are the total pressure and the total temperature 
at the block inflow face. In addition, this particular 
inflow boundary condition assumes a fluid flow an- 
gle normal to the inflow face. The operating NPR 
of the nozzle and the frec-stream static pressure p x 
determine the jet total pressure at the inflow face 
from p tJ = (NPR)(p x ). The design NPR is calcu- 
lated from the series of equations shown in appon- 
dix A. A typical jet total temperature obtained from 
the experimental tests is used for all the calculat ions. 

A characteristic boundary condition with 
Riemann invariants is applied to external inflow and 
outer lateral boundaries. An extrapolation bound- 
ary condition is applied on the downstream outflow 
face where the nozzle plume exits the computational 
domain. 

Performance Calculation 

Nozzle performance is obtained through the ap- 
plication of the momentum theorem to a control 
volume surrounding the nozzle (ref. 2). Cheatham, 
Walker, and Gridley calculated both 2-D and ,'M) in- 
viscid nozzle performance using this method (ref. 9). 
The surface over which the integration of the flow 
quantities is performed is typically the nozzle exit. 
The method utilized herein integrates the mass and 
momentum fluxes and the pressure forces over each 
incremental cell modeling the nozzle outflow faces 
through the use of equations (1) and (2): 

W P = £ (/>U • n ) A4 (1) 

F = E [^ U ( u ' n ) + (p - Ptc) n] A A (2a) 

M = F x R (2b) 
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where A A is the outflow area attributed to the 
incremental cell face and R is the moment arm from 
the reference center to the incremental cell. 

Ideal mass-flow rate and thrust are determined 
from the isentropic flow equations (eqs. (3) and (4)) 
and are used to normalize the calculated mass-flow 
rate and thrust for comparisons with the experimen- 
tal data: 
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The performance package is incorporated into 
PAB3D to permit monitoring of various performance 
parameters as the solution converges. 

Results and Discussion 

Solutions were obtained for four nozzle configu- 
rations from the flow method PAB3D. The solutions 
were all computed with laminar viscous stress mod- 
eling due to biasing of the turbulent flow models 
to a single direction. This modeling obviates tur- 
bulence calculations when an H-type grid is used to 
describe the nozzle geometry. Discharge coefficient, 
thrust ratio, and off-axis force and moment quanti- 
ties (where applicable) were calculated for different 
operating conditions. 

Case Is Stratford Choke Nozzle Study 

A sectional view of the Stratford choke nozzle is 
shown in figure 2. Gridding arrangements and con- 
vergence data arc shown in figures 7 to 10. Four com- 
putational grids were run for this geometry. Residual 
trends and discharge coefficient convergence are plot- 
ted for each grid. 

Grid definition. This case study was modeled 
with two computational zones. The first zone con- 
tained the nozzle internal grid and an external grid 
region surrounding the nozzle, and the second zone 
modeled the external air downstream of the noz- 
zle exit. The dimensions of the nozzle block were 
129 x 2 x 105. The nozzle internal-boundary-layer 
grid region contained 34 points, with the first grid 
point specified from = 4. Grid expansion rates 
were less than 5 percent in the axial and radial di- 
rections. The external region surrounding the noz- 
zle grid was dimensioned at 129 x 2 x 208. Total- 
pressure and total-temperature boundary conditions 


were applied to both inflow faces of these blocks. 
The external downstream region was dimensioned at 
145 x 2 x 313. The zone interface matched point for 
point between the two zones. The grid was stretched 
toward the outer boundary to a distance of 40 exit 
radii. Solid-wall boundary conditions were applied 
to this surface. The outflow boundary was 15 exit 
radii downstream of the nozzle exit plane. Extrapo- 
lation boundary conditions were applied to the out- 
flow boundary. Each of the successively coarser grids 
used grid points from this mesh. A sectional view 
of the grid in the zfc-plane is shown in figure 7(a). 
This grid represented the most dense concentration 
of points for this study. The grid density was reduced 
by two in both the i-direction (axial), and /c-direction 
(radial) three successive times. Figures 8(a), 9(a), 
and 10(a) are sectional views of grids reduced by V2, 
1 / 4 , and Vs. Grid 1 is the baseline grid. Grids 2, 3, 
and 4 are the V2, V4, and Vs grids, respectively. 

An ax isym metric flow symmetry is assumed 
through the use of a 2-D wedge-angle boundary con- 
dition. A single-cell-width polar grid was generated 
with a wedge angle of 5.625°, which represents Vg 4 of 
an axisymmetric geometry. The wedge-angle bound- 
ary condition calculates the angle of the wedge lateral 
faces and determines the proper boundary values to 
apply to satisfy the symmetry assumptions. 

Discharge coefficient correlation. The noz- 
zle pressure ratio for this study was set at 2. I he 
external free stream was initially set at A/x — 0.05. 
The initial solutions for the time-dependent calcula- 
tions for grids 2, 3, and 4 were produced through a 
two-step procedure. A one-dimensional (1-D) isen- 
tropic solution was calculated for the nozzle internal 
grid, as described in appendix B. This flow field was 
the starting point for a space- marching solution for 
the jet exhaust downstream of the nozzle exit. Time- 
dependent calculations were continued on the result- 
ing flow field. The calculation for grid 1 utilized a 
converged solution from grid 2 extrapolated to the 
grid 1 computational domain for the starting solu- 
tion. Residual trends for the four grids are shown in 
figures 7(b), 8(b), 9(b), and 10(b). Discharge coef- 
ficient convergences are shown in figures 7(c), 8(c), 
9(c), and 10(c). The effect of grid count on discharge 
coefficient is shown figure 11. 

The discharge coefficient stabilized for grids 2 
and 3 within 1000 solution iterations (1 sweep down- 
stream and 1 sweep back upstream by the solver 
was counted as 1 iteration) and grid 4 stabilized 
within 600 iterations. There appears to have been 
no close correlation between the solution residuals 
and convergence of the calculated nozzle discharge 
coefficient after the first 1000 iterations. Excursions 
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of 2 orders of magnitude in the residual were ob- 
served, with no discernible change in the discharge 
coefficient after that point. The discharge coefficient 
for grid 1 remained stable as well, despite the resid- 
ual dropping over an order of magnitude over ap- 
proximately 1000 iterations. Solution residuals were 
therefore used only as an initial indicator of solution 
stability (i.e., that the solution is not diverging) and 
discharge coefficient was utilized as the indicator of 
final solution convergence. 

The effect of grid count of discharge coefficient is 
shown in figure 11. The very coarse grid (grid 4) of 
this study converged to 0.5 percent of the level of 
the experimental data. Grid 3 converged to within 
0.1 percent, and grids 2 and 1 were less than 0.1 per- 
cent of the experimental discharge coefficient. Dis- 
charge coefficient was a stable and convergent quan- 
tity with grid count. 

Case 2: Throat Contouring Study 

Grid definition . Two nozzles (Al and A2, see 
fig. 3) with different throat radii were modeled. A 
representative internal computational grid for a non- 
axisymmetric nozzle is shown in figure 12(a). The 
ordinate of the plot was nondirnensionalized by the 
length of the divergent flap, with the origin at the 
nozzle throat. The computational grid for the nozzle 
was an H-H grid consisting of 61 grid points in the 
^-direction (axial), 49 points in the j-direction (lat- 
eral), and 37 grid points in the A>direction (vertical). 
The grid spacing of the internal grid mesh was packed 
in the boundary layer near the walls and stretched 
away in the center region of the nozzle. Flow sym- 
metry was assumed about the horizontal and vertical 
planes to reduce computational time. The number of 
cells was identical for both nozzle configurations. 

The nozzle internal-boundary-layer grid region 
contained 13 grid points, with the first grid point 
specified as 0.00003 cm from the wall surface. This 
region was roughly 0.115 cm deep in the fc-direction 
along the divergent flaps and 0.251 cm deep in the 
j-direction along the sidewalls at the nozzle exit 
(exit half-height was roughly 1.5 cm). Constant- 
interval spacing from the boundary-layer edge to the 
centerline was used for the remaining grid points. 

In the calculation of off-design settings (i.e., 
underexpanded nozzle flow), the static external air 
was modeled to allow the jet plume to adjust to 
the jump in static pressure that occurred near the 
nozzle exit face. The external and internal domains 
were split into three zones, as shown in figures 12(b) 
and 12(c). The first zone (61 x 49 x 37) contained 
the nozzle internal geometry, which was the same as 


in the single-zone case. The second (8 x 87 x 73) 
and third (17 x 44 x 37) zones were free-stream do- 
mains, with the second zone serving as a transition 
between the denser nozzle geometry grid and the rel- 
atively sparse grid of the third zone. This transition 
avoided radical grid changes, such as the 2:1 grid 
reduction at the boundary to the second and third 
zones, in areas of possibly high flow gradients such 
as those expected to occur at the nozzle exit. A rep- 
resentative grid arrangement for the first z-plane of 
the second zone is shown in figure 12(d). A cut in 
this zone face matched point for point the nozzle exit 
grid of the last f-plane in the first zone. The grid was 
stretched toward the outer boundaries at a distance 
of five equivalent nozzle radii away from the outer 
line of the matching cut. The code required about 
4 million words of memory for the execution of the 
multizone cases (approximately 190 000 grid points). 

A 1-D isentropic solution described in appendix B 
was used for the starting flow field. Mass flow of the 
nozzle was observed to gauge solution convergence. 

Pressure data correlation . The computed 
nozzle internal pressures for nozzles Al and A2. 
which arc presented in figure 13, compared well with 
experimental data. These solutions converged within 
about 400 time steps that took roughly 3100 sec 
on the Cray-2 computer with an initial Courant- 
Friedrichs-Lewy (CFL) stability number of 3. The 
CFL numbers for subsequent iterations during a 
run are adjusted automatically within the code to 
accelerate convergence. 

The ability to calculate the change in nozzle flow 
due to a difference in throat radius (static pressure 
distribution of nozzle Al subtracted from that of 
nozzle A2) is shown in figure 14. The magnitude and 
location of the incremental pressure differences due 
to the change in nozzle throat radius were fairly well 
matched. This pressure difference integrated over 
the length of the divergent flaps (i.e., x/l s > 0) was 
1.58 N, or roughly an increment of 0.002 in thrust 
ratio based on ideal thrust (from the ideal mass flow). 
The pressure gradients across the exit were lower for 
nozzle A2 than for nozzle Al, probably because of the 
weaker pressure gradients between the divergent flaps 
of nozzle A2. The magnitude of the pressure force 
at the exit (subtracted from the free-stream static 
pressure) was much smaller for nozzle A2 than for 
Al (-0.47 N compared with -11.12 N). It appears 
that nozzle A2 was operating closer to the theoretical 
design conditions than nozzle Al was operating for 
the same input total-pressure conditions because of 
the smaller loss of thrust to pressure forces. 
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Nozzle performance . The discharge coefficient 
converged to within 0.5 percent of the experimental 
data level. This level of convergence was typical for 
this particular set of single-zone calculations of the 
nozzle internal flow. 

As shown in figure 15, the agreement of the com- 
puted and the experimental discharge coefficient was 
within experimental accuracy for both nozzles. Fig- 
ure 16 shows the comparison of calculated and mea- 
sured internal thrust ratios for the range of NPR’s 
tested for nozzles Al and A2. The calculated oil- 
design (NPR r y = 2.97) thrust ratio was 0.6 percent 
high for nozzle Al and IT percent high for nozzle A2 
as compared with the experimental data. This com- 
parison was relatively consistent over the range of 
NPR’s from design to 8.0. Laminar viscous stress 
modeling could have contributed to excess thrust be- 
ing calculated because a boundary layer that was 
thinner than actual was calculated along the diver- 
gent flaps. 

With a Reynolds number based on the equiva- 
lent throat radius (approximately 1.28 x 10 f> ). lam- 
inar (from Blasius’ equation) and turbulent (from 
ref. 10. p. 636) boundary- layer profiles were calcu- 
lated for flow over a flat plate for a distance' equal 
to the length of the nozzle divergent flaps. A por- 
tion of these profiles and the velocity distribution 
at the exit plane calculated by PAB3D for noz- 
zle Al arc' shown in figure 17. As would be expected 
from a flow calculat ion using laminar viscous st resses 
in PAB3D, the laminar-flat-plate distribution more 
closely matched the boundary layer calculated by 
PAB3D than did the turbulent-Hat-plate profile'. The 
laminar-boundary-layer thickness occupied roughly 
2 percent of the nozzle semi-height at the exit. A sim- 
ilar calculation using a turbulent flow model resulted 
in a boundary-layer thickness of almost 9 percent of 
the nozzle exit semi-height. The velocity distribu- 
tion along the upper divergent flap showed several 
regions of reversed flow, as shown in figure 18. The 
first boundary-layer reattach ment, a few centimeters 
downstream of the throat, was a likely location of 
the transition to turbulent How. Therefore, the ac- 
tual boundary-layer thickness at the nozzle exit may 
have been slightly thinner than that of the turbu- 
lent flow estimate, but thicker than that of just the 
laminar calculation. A thicker boundary layer on the 
nozzle divergent flaps would result in lower computed 
values of thrust. 

The difference in performance between nozzles Al 
and A2 is shown in figure 19 for experimental and 
theoretical results. The difference in experimen- 
tal nozzle thrust performance due to the change in 
throat radius was statistically zero, as the variance in 


performance was less than the experimental data ac- 
curacy of 0.5 percent. Similar small increments were 
computed by PAB3D for the range of NPR’s from 
design to 8.0. 

Case 3: Translating Sidewall Yaw 

Vectoring Study 

A cross-sectional sketch of the translating sidewall 
yaw vectoring concept of reference 5 is shown in fig- 
ure 4. Effective yaw vector angle was experimentally 
measured for five sidewall positions at several NPR 
settings. The performance validation for this set of 
nozzle geometries cut across the experimental data 
two ways. The first fixed the geometry at one side- 
wall position, with calculations made for three pres- 
sure ratios: design NPR (NPR,/ = 2.967), and NPR 
— 5.0 and 8.0. The second was made on geometries 
of varying sidewall cutback positions, with the nozzle 
operating at a constant NPR = 5.0. The quarter- 
plane symmetric-flow model could no longer be used 
because of the generation of unbalanced off-axis (i.e., 
yawing) forces. The generation of side forces with 
no net normal forces would have required this sym- 
metry boundary to be placed in the horizontal plane 
cutting through the center of the model. 

Grid definition . The computational grid was 
split into 4 streamwise zones, with the first zone con- 
taining 2 blocks and the remaining 3 zones containing 
10 blocks each. A representative sketch of the grid is 
shown in figure 20. A quarter-plane cross-sectional 
cut of the grid at the nozzle exit plant? and the block- 
ing strategy art' shown in figures 20(c) and 20(d). 
Blocks 1 and 2 in zone 1 modeled the internal flow 
path containing 41 x 49x37 grid points, each using an 
H-H topology. The cross-plane and boundary-layer 
grid was similar to that used for case 2. More densely 
packet 1 grid is placed in the external lateral sidewall 
region, where the vectored flow could exhaust into 
tht? static free stream. All the external blocks used 
H-O grid topology. 

Performance calculation . Calculation of ef- 
fective yaw vector angle required the determination 
of both the axial force and the side force generated 
by the nozzle flow. Integration of the momentum 
fluxes and the shearing forces across the nozzle exit 
face and the cutaway part of the sidewall accounted 
for all the vectoring forces. The prediction of the 
change in effective yaw vector angle with NPR for a 
fixed sidewall cutback position is shown in figure 21. 
Prediction of the change in effective yaw vector an- 
gle with sidewall position for a fixed NPR is shown 
in figure 22. 
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Overall, the trend and the level of effective yaw 
vector angle were predicted to within 0.5° of the 
experimental data over the range of NPR’s for noz- 
zle S3 ((x. s — Xf )/l s = 0.25). The static pressure at 
the nozzle exit is lower than the free-stream static 
pressure when operating at overexpanded conditions. 
This pressure difference caused air to be drawn into 
the nozzle passage through the cutback sidewall; 
hence, the sign of the effective yaw vector angle; at 
the low NPR was reversed. Operating the nozzle at 
design conditions resulted in a slightly positive yaw 
vector angle because of a small amount of flow out- 
ward through the sidewall cutout. The zero point for 
the effective yaw vector angle would have occurred at 
some operating condition between design NPR and 
an overexpanded NPR, where the momentum fluxes 
creating side forces balanced out across the available 
outflow surfaces of the nozzle. 

The change in effective yaw vector angle with 
sidewall cutback was predicted to within 0.3° over 
the range of sidewall positions at NPR = 5.0. (See 
fig. 22.) Along with the integration of the mo- 
menta flux and pressure forces in the lateral direc- 
tion through the sidewall cutout, the shearing forces 
across the exit face must be included to account for 
all the sideward forces. The shearing forces for noz- 
zle S3 contributed 75 percent of the total side force 
generated by the nozzle. Similarly, the axially di- 
rected shearing forces along the open sidewall face 
contributed roughly 10 percent to the thrust of the 
nozzle. 

Case 4: Multiaxis Thrust Vectoring Study 

The axisymrnetric nozzle geometries reported by 
Carson and Capone (ref. 6) consisted of designs for 
two power settings and two flap lengths. (See fig. 5.) 
The short-flap, dry power nozzle with NPR,/ = 5.01 
(e = 1.350) was chosen for this set of comparisons. 
This nozzle was tested at four geometric pitch vec- 
tor angles (# t; , p = 0°, 10°, 20°, and 30°) and one 
combined pitch-yaw vectoring mode (S v , p = 14° , 
fiv,y = 14°). A sectional sketch of the nozzle internal 
shape for the different pitch vector angles is shown in 
figure 6. Yawing and pitching moments along with 
mass flow and axial thrust were experimentally mea- 
sured. The moments were nondimensionalized by the 
ideal axial thrust and a reference length. The refer- 
ence length was the nozzle throat diameter, 5.756 cm, 
and the moment reference center was the center of 
the force balance, 24.765 cm upstream of the nozzle 
connect station. Predicted and experimental data 
are presented for the nozzle operating at the design 
NPR. 


Grid definition . The prediction of the off-axis 
forces required a gridding of the geometry with no as- 
sumed planes of symmetry. The computational do- 
main was divided into two zones, with two blocks 
in each zone for the nozzle split laterally. An oblique 
view of the cell-centered grid is shown in figure 23(a). 
The grid in the centerline plane is shown in fig- 
ure 23(b) for the nozzle vectored 0° and in fig- 
ure 23(c) for the nozzle vectored 30°. The first zone, 
with each block dimensioned 28 x 17 x 17, extended 
from just upstream of the physical nozzle connect 
station to the beginning of the convergent section 
before the throat. The second zone, with each block 
dimensioned 45 x 33 x 33, consisted of the conver- 
gent and divergent sections of the nozzle ending at 
the nozzle exit plane. Each block had an H-O grid 
topology, as shown in figures 23(c) and 23(d). An 
extrapolation boundary condition applied to the exit 
face allowed for terminating the computational vol- 
ume at that point for calculations of flow at on-design 
nozzle pressure ratios. 

Performance calculation . Experimental data 
and various calculated performance quantities for 
four different configurations are shown in figures 24 
to 27. The numerical schemes used for solutions for 
this geometry were first order in the ^-direction and 
third order in the i- and indirections. All five con- 
figurations were run using laminar flow assumptions. 
Nozzle mass flow was determined from the statistical 
average of the mass flow through the first 18 /-planes 
of the first zone. The moments and forces generated 
by the nozzle flow were obtained through the inte- 
gration of the fluid momentum flux, pressure, and 
shearing forces at the exit plane. 

The level of and the increment in discharge coeffi- 
cient with change in the geometric pitch vector angle 
were predicted to within experimental accuracy at 
the design nozzle pressure ratio of 5.01. (See fig. 24.) 
Losses due to the sharp throat contour kept w p /wj at 
around 0.95 for f) v . p = 0° and 10°. Gradually larger 
losses were predicted for the higher vectoring angles. 

The change in thrust ratio with varying geometric 
pitch angle was predicted to within 2 percent of 
the experimental data at, the design nozzle pressure 
ratio, as shown in figure 25. The change in pitching 
moment with varying geometric pitch vector angle 
was similarly predicted, as shown in figure 26. 

The capability of predicting multiaxis thrust vec- 
toring is shown in figure 27 for the following three 
configurations: S 1Kp = 0°, 6 vpj = 0°; 6 v:p = 20°, 
6 v ,y ~ 0°; and b vp — 14°, b luy = 14°. Experi- 
mentally, the combined pitch-yaw thrust vectoring 
case was generated by a 45° rotation about the 
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centerline axis at the nozzle connect station of the 
nozzle with 20° of pitch vectoring. Exactly the same 
methodology was numerically applied to the ana- 
lytical description of the 20° pitch vector geometry 
to produce the combined pitch-yaw computational 
grid. Numerically the yawing moments predicted 
for the 6 v ^y = 0° geometries were zero. The pre- 
dicted yawing moment for the combined pitch-yaw 
geometry was within 0.7 percent of the experimental 
data (fig. 27(b)). Similar agreement of experimen- 
tal data with predicted pitching moment is shown in 
figure 27(a). The relative levels of pitching moment 
were within 1 percent of experimental levels, and in 
particular the loss of pitching moment due to rota- 
tion of the pitch axis by 45° was predicted. 

Conclusions 

A nozzle internal performance module has been 
written and incorporated into a three-dimensional 
Navier-Stokes solver PAB3D for the calculation of 
nozzle internal performance. The flow quantities 
calculated are discharge coefficient; normal, side, 
and axial thrust ratios; rolling, pitching, and yawing 
moments; and effective pitch and yaw vector angles. 
The calculations were compared with experimental 
data from several investigations performed in the 
static test facility of the Langley 16-Foot Transonic 
Tunnel. The results are summarized as follows: 


1. Nozzle discharge coefficient was predicted to 
within experimental accuracy over a range of noz- 
zle pressure ratios from design to underexpanded 
pressure ratios. 

2. The effect of throat contour geometry on noz- 
zle discharge coefficient and thrust ratio was pre- 
dicted over a range of nozzle pressure ratios for a 
two-dimensional convergent- divergent nozzle. 

3. The effect of sidewall cutback on the effective 
yaw vector angle generated by a two-dimensional 
convergent-divergent nozzle operating at an 
underexpanded pressure ratio was predicted typi- 
cally to within 0.5° of the experimental data. The 
effect of varying nozzle pressure ratio from over- 
expanded to underexpanded conditions on effec- 
tive yaw vector angle for a fixed sidewall cutback 
position was also predicted typically to within 
0.5 percent of the experimental data. 

4. Thrust vectoring forces and moments were pre- 
dicted typically within 2 percent of experimental 
data for an axisymmetric multiaxis thrust vector- 
ing nozzle operating at the design nozzle pressure 
ratio. 

NASA Langley Research Center 

Hampton, VA 23681-0001 

August 13, 1992 
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Appendix A 

Calculation of Design Nozzle Pressure 
Ratio 

The physical geometry of a nozzle fixes the noz- 
zle pressure ratio (NPR) at which the most efficient 
operation of the nozzle occurs. The following equa- 
tions are used to determine the design NPR when the 
nozzle expansion ratio e is known. 


by using a Newton iterative method on the following 
equation: 



(A2) 

where ^ is used for The ratio T j s calculated 
from the Mach number: 


The inverse of the nozzle expansion ratio is the 
ratio of the nozzle throat area to the nozzle exit area: 

A t 1 

A: = 7 < A1 > 

The nozzle exit Mach number (isentropic) is solved 


P_ 

Pt 




(A3a) 


The inverse of 
nozzle: 


R 

Pt 


is the design pressure ratio of the 

NPR '< - Jz <A3b) 
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Appendix B 

Calculation of Initial Flow-Field 
Solutions 

A starting Mach number distribution is deter- 
mined from the area distribution of the internal ge- 
ometry using equation (A2) and assuming sonic 

conditions at the throat. The dimensionless ratios of 
temperature and pressure are then calculated from 
the Mach number distribution interpolated to cell- 
centered coordinates with the following 1-D isen- 
t.ropie relations: 


T 

T, 



(Bl) 


= (^ + 1— U/, 2 V ' (B2) 

pt V 2 / 

Typical scale-model operating jet total-pressure 
and total-temperature levels are chosen to calculate 


the quantities p, 
tions: 


u, and e with the following equa- 


II 

c.. 

(B3) 

P 

P = PtJ ~ 

Pi 

(B4) 

P 

P RT 

(B5) 

a = y/fRT 

(BO) 

a = ( iM i 

(B7) 

1 2 , V 

= 2 PU + 7 - 1 

(B8) 


Eacli plane of c*t'lls (/ = Constant ) is uniformly as- 
signed these quantities. The initial solution assumes 
the cross-components of velocity (e.g.. v and w) are 
zero. The flow variables p, pu, pv , pw, and e are then 
written to the restart file of the code and utilized as 
the initial solution. 
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Residual 



(b) Variation of solution residual with solution sweep number 


Cd 


1.0300 

1.0200 

1.0100 

1.0000 

.9900 

.9800 

.9700 

.9600 

.9500 

.9400 

.9300 


500 1000 1500 2000 2500 3000 

Sweep number 


3500 4000 4500 


(c) Convergence of discharge coefficient with solution sweep number. 

Figure 7. Concluded. 
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Nozzle exit 


(a) Grid showing ifc-plane. 

Figure 8. Grid density reduction of 'd for Stratford choke nozzle (case 1). 
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Residual 
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(b) Variation of solution residual with solution sweep number. 



(c) Convergence of discharge coefficient with solution sweep number. 

Figure 8. Concluded. 
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Nozzle exit 

(a) Grid showing ik- plant'. 

Figure 9. Grid density reduction of >A for Stratford choke nozzle (case 1). 
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Residual 



(b) Variation of solution residual with solution sweep number. 



(c) Convergence of discharge coefficient with solution sweep number. 

Figure 9. Concluded. 
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Residual 
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Sweep number 

(b) Variation of solution residual with solution sweep number. 



(c) Convergence of discharge coefficient with solution sweep number. 
Figure 10. Concluded. 
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□ Upper flap (ref. 4) 
O Lower flap (ref. 4) 
— Theory, PAB3D 



(a) Nozzle configuration Al. 



(b) Nozzle configuration A 2. 

Figure 13. Static pressure distributions of theory and experiment, on upper-flap centerline for 
power; design NPR. 


case 2. Dry 


.200 


O Experiment (ref. 4) 
Theory, PAB3D 



Figure 14. Experimental and theoretical increments in static pressure on upper-flap centerline due to change 
in nozzle throat radius for case 2. Dry power; design NPR. 



NPR 


Figure 15. Predicted and experimentally determined discharge coefficients for case 2. 
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O Experiment, nozzle A1 (ref. 4) 
Theory, PAB3D 



O Experiment, nozzle A2 (ref. 4) 
Theory, PAB3D 



NPR 

(b) Nozzle A2. 

Figure 16. Experimental and theoretical internal thrust ratios for case 2 




Figure 17. Calculated boundary-layer profile on divergent flap at nozzle exit for nozzle configuration A1 









Nozzle exit 
plane 


(b) Vertical sectional view of centerline plane of symmetry. 
Figure 20. Continued. 
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(c) Vertical sectional view in plane of nozzle exit, one-quarter grid in lateral ( ZY ) plane. 

Figure 20. Continued. 
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NPR 

Figure 21. Predicted and experimentally determined changes in effective yaw vector angle with NPR for case 3. 
Dry power nozzle; (x <s — xt)/l H = 0.25. 



Figure 22. Predicted and experimentally determined changes in effective yaw vector angle with sidewall position 
for case 3. Dry power nozzle; NPR = 5.0. 



(b) Vertical centerline plane of symmetry of unvectored nozzle for short-flap, dry power geometry. 



(c) Vertical centerline plane of symmetry of nozzle vectored 30° for short-flap, dry power geometry. 

Figure 23. Continued. 
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(d) Cross-sectional view in plane of nozzle exit face. 
Figure 23. Concluded. 
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(b) Yawing moment. 


Figure 27. Concluded. 
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